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Abstract 

The 6-hadron production cross section is measured with the ATLAS detector in pp colli- 
sions at y/s = 7 TeV, using 3.3 pb -1 of integrated luminosity, collected during the 2010 LHC 
run. The 6-hadrons are selected by partially reconstructing D* + fi~X final states. Differential 
cross sections are measured as functions of the transverse momentum and pseudorapidity. 
The measured production cross section for a 6-hadron with p^ > 9 GeV and \r]\ < 2.5 is 
32.7 ± 0.8 (stat.)^gg (syst.) fib, higher than the next-to-leading-order QCD predictions but 
consistent within the experimental and theoretical uncertainties. 
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1. Introduction 

The production of heavy quarks at hadron colliders provides a challenging opportunity 
to test the validity of quantum chromodynamics (QCD) predictions and calculations. The 
6-hadron production cross section has been predicted with next-to-leading-order (NLO) ac- 
curacy for more than twenty years [1, 2j. 

Several measurements were performed with proton-antiproton collisions by the UAl ex- 
periment at the SppS collider (CERN) at a centre-of-mass energy of s/s = 630 GeV [31 Hj, 
and by the CDF and DO experiments at the Tevatron collider (Fermilab) at a/s = 630 GeV, 
1.8 TeV and 1.96 TeV [5T-TH]. These measurements made a significant contribution to the 
understanding of heavy-quark production in hadronic collisions [15] , but the theoretical 
predictions still suffer from large uncertainties, mainly due to the dependence on the fac- 
torisation and renormalisation scales. 

A measurement of the 6-hadron production cross section in proton-proton collisions at the 
Large Hadron Collider (LHC) provides a further test of QCD calculations for heavy-quark 
production at higher centre-of-mass energies. Recently the LHCb experiment measured the 
bb and B + [TBTTIS] production cross sections in the forward region at -y/s = 7 TeV, the 
CMS experiment measured the production cross sections for B + , B°, B° s mesons, inclusive 
6-hadrons with muons, and bb decays with muons at a/s = 7 TeV [T9H23] . and the ALICE 
experiment measured the bb production cross section in pp collisions at a/s = 7 TeV [21] . 

This paper presents a measurement of the 6-hadron (H^, a hadron containing a 6-quark 
and not a 6-quark) production cross section at a centre-of-mass energy of 7 TeV with the 
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ATLAS detector at the LHC, and its comparison with the NLO QCD theoretical pre- 
dictions. The measurement requires the partial reconstruction of the 6-hadron decay fi- 
nal state D* + fi~X, with the D* + reconstructed through the fully hadronic decay chain 
D* + — > n + D°(— K~n + ). This sample was collected by ATLAS between August and Octo- 
ber 2010 using events selected by a single-muon trigger, and corresponds to a total integrated 
luminosity of 3.3 pb -1 . 

2. The ATLAS detector 

The ATLAS detector (25] covers almost the full solid angle around the collision point 
with layers of tracking detectors, calorimeters and muon chambers. For the measurement 
presented in this paper, the inner detector tracking devices, the muon spectrometer and the 
trigger system are of particular importance. 

The inner detector (ID) has full coverage in <fr and covers the pseudorapidity range 
\r]\ < 2.5. It consists of a silicon pixel detector, a silicon microstrip tracker and a transition 
radiation tracker composed of drift tubes. These detectors are located at radial distances 
of 50.5-1066 mm from the interaction point and are surrounded by a thin superconducting 
solenoid providing a 2 T axial magnetic field. The ID barrel consists of three layers of pixels, 
four double-layers of single-sided silicon microstrips, and 73 layers of drift tubes, while each 
ID end-cap has three layers of pixels, nine double-layers of single-sided silicon microstrips, 
and 160 layers of drift tubes. 

The muon spectrometer covers the pseudorapidity range \r]\ < 2.7 and is located within 
the magnetic field produced by three large superconducting air-core toroid systems. The 
muon spectrometer is divided into a barrel region < 1.05) and two end-cap regions 
(1.05 < \r]\ < 2.7), within which the average magnetic fields are 0.5 T and IT respectively. 
Precise measurements are made in the bending plane by monitored drift tube chambers, 
or, in the innermost layer for 2.0 < < 2.7, by cathode strip chambers. Resistive plate 
chambers in the barrel and thin gap chambers at < 2.4 in the end-caps are used as trigger 
chambers. The chambers are arranged in three layers, such that high muons traverse at 
least three stations with a lever arm of several metres. 

A three-level trigger system is used to select interesting events. The first level is hardware- 
based, and uses a subset of the detector information to reduce the event rate to a design 
value of at most 75 kHz. This is followed by two software-based trigger levels, together 
known as the high level trigger, which finally reduce the event rate to about 200 Hz. 

3. Outline of the measurement 

The first result presented in this paper is the Hi, — > D* + X production cross section, 
measured in a limited fiducial acceptance for the D* + /i~ final state. Given the integrated 
luminosity J?f of the data sample, and the branching ratio SS of the D* + cascade decay 
D* + — > n + D°(— > K~n + ), the H b — > D* + ^rX cross section is defined as: 
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where N(D* + yT + D*~ is the total number of reconstructed candidates, /& is the fraction 
of candidates originating from the decay Hj, — > D* + fi~X and e is the signal reconstruction 
efficiency. The efficiency takes into account reconstruction and muon trigger efficiencies, 
including the loss of events where the D* + falls within the fiducial acceptance, but the 
decay products (it or K) can not be reconstructed because they fall outside the pt and 77 
acceptance. The number N of reconstructed candidates includes both D* + pT and D*~ ji + 
combinations: assuming that b- and 6-quarks are produced with the same rate at the LHC, 
the factor of two is needed to quote the cross section for hadrons containing a 6-quark. The 
value of the branching ratio BB can be obtained by combining the world average values of 
the branching ratios D* + -»■ n+D° and D° K~ir + [26], and is (2.63 ± 0.04) %. 

The parameters N, fa and e are determined as functions of the transverse momentum 
and pseudorapidity of the D* + fi~ pairs, in order to measure the differential cross sections. 
The detailed calculation of these parameters is discussed in the following sections. 

To obtain the fe-hadron production cross section a(pp —> HbX), the — > D* + n~ X 
cross section is divided by an acceptance correction a, accounting for the fiducial region in 
which this is measured, and by the inclusive branching ratio BB(b — > D fi~X). For this 
branching ratio the world average value is (2. 75 ±0.19) %, assuming the world average values 
of the 6-hadronisation fractions [26]. The dominant contributions to the sample are from 
B° mesons, through the decay B° — > D*~ /i + z/ M and its charge conjugate. 

4. Event simulation and NLO cross section predictions 

Monte Carlo (MC) simulated samples are used to optimise the selection criteria (Sec- 
tion p| and to evaluate the D* + fi~ signal composition and reconstruction efficiency (Sec- 
tions^6] and [7|) . The different b- and c-quark sources of D* + \T are studied using inclusive 
samples of bb and cc events having at least one muon with > 4 GeV and \r]\ < 2.5 in 
the final state. Both samples are generated with Pythia [27], using the ATLAS AMBT1 
tuning [28]. The ATLAS detector response to the passage of the generated particles is simu- 
lated with Geant4 [221130], and the simulated events are fully reconstructed with the same 
software used to process the collision data. 

To compare the measurements with theoretical predictions, NLO QCD calculations, 
matched with a leading-logarithmic parton shower MC simulation, are used. Predictions 
for bb production at the LHC at a/s = 7 TeV are evaluated with two packages: mc@nlo 
4.0 [2DE2] and Powheg-hvq 1.01 [331EII- mc@nlo is matched with the Herwig 6.5 [33] 
MC event generator, while Powheg is used with both Herwig 6.5 and Pythia 6.4 [27] . 
For all the predictions, the inclusive branching ratio SSib — > D* + fi~X) is set to the world 
average value. 

The following set of input parameters is used to perform all theoretical predictions: 

• CTEQ6.6 |36j parameterisation for the proton parton distribution function (PDF). 

• 6-quark mass mj of 4.75 GeV |26j . 
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• Renormalisation and factorisation scales set to /x r = /iy = fi, where /x has different 
definitions for MC@NLO and POWHEG. For MC@NLO: 



2 2 i 

fi = m Q + 



(Pt,q+Pt,q) 2 



where Pt,q arid Pt,q are the transverse momenta of the produced heavy quark and 
antiquark, and mg is the heavy-quark mass. For Powheg: 

^ = m Q + ( m Qg/ 4 - m o) s ^( 6 q) 

where mgq is the invariant mass of the QQ system and 9q is the polar angle of the 
heavy quark in the QQ rest frame. 

• Heavy-quark hadronisation: cluster model [37J for Herwig; Lund string model [38] 
with Bowler modification [39] of the Lund symmetric fragmentation function [3D] for 
Pythia. 

The following sources of theoretical uncertainties are included in the NLO predictions: 

• Scale uncertainty, determined by varying \i r and /i/ independently to ji/2 and 2/x, 
with the additional constraint 1/2 < /i r //J/ < 2, and selecting the largest positive and 
negative variations. 

• nib uncertainty, determined by varying the 6-quark mass by ±0.25 GeV. 

• PDF uncertainty, determined by using the CTEQ6.6 PDF error eigenvectors; the total 
uncertainty is obtained by varying each parameter independently within these errors 
and summing the resulting variations in quadrature. 

• Hadronisation uncertainty, determined in Pythia by using the Peterson fragmen- 
tation function [41] instead of the Bowler one, with extreme choices of the 6-quark 
fragmentation parameter: e& = 0.002 and ej, = 0.01. 

In addition to the final comparison with the experimental measurement, these theoretical 
predictions are used to unfold and extrapolate the measured cross sections (Sections M and 



10), and to extrapolate to the full kinematic phase space (Section [TT]). In the following 



Powheg+Pythia is used as the default prediction. 



5. Data selection and reconstruction of the D* + fi decay 

The D* + ur (including its charge conjugate) sample was collected during stable proton- 
proton collisions. Events were selected by a single-muon trigger, which requires a muon, 
reconstructed by the high level trigger, with p-? > 6 GeV. This trigger was prescaled during 
the last part of the 2010 data-taking period. Taking into account the prescale factors, this 
data sample corresponds to an integrated luminosity of 3.3 pb -1 . 
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The D* candidates are reconstructed through the fully hadronic decay chain D* — > 
n + D (— > K~n + ), using only good quality tracks, i.e. tracks with at least five silicon detector 
hits, and at least one of them in the pixel detector. 

The 6-hadron and D° decay vertices are reconstructed and fitted simultaneously. To 
perform the vertexing, an iterative procedure based on a fast Kalman filtering method is 
used. This allows to reconstruct consecutively all the vertices of the same decay chain, 
using the full information from track reconstruction (particles trajectories with complete 
error matrices). All pairs of opposite charge particle tracks are fitted to a single vertex to 
form D° candidates, assigning to each track, in turn, the kaon or the pion mass, with the 
additional requirement px > 1 GeV for both the kaon and pion candidate; the resulting D° 
candidate is reconstructed by combining the kaon and pion four-momenta. The D° path is 
then extrapolated back and fitted with a track of opposite charge to the candidate kaon, 
requiring px > 250 MeV and assigning to it the pion mass, to form the D* + candidate, and 
with a muon with px > 6 GeV and \q\ < 2.4 to form the 6-hadron vertex. No requirements 
are made here on the muon charge; only opposite charge combinations D* + fi~ are used 
in the analysis, while same charge combinations are used to cross-check the background. 
The muon is also required to have fired the trigger. To ensure good fit quality, the global 
X 2 probability of the combined fit must satisfy P(x 2 ) > 0.001. To avoid an additional 
systematic uncertainty no requirement on the 6-hadron vertex decay length is applied. 

The D* + candidate is accepted if it satisfies Pt{K ~7r + 7r + ) > 4.5 GeV and |^(K _ 7r + 7r + ) | < 
2.5, and either (a) \m(K~n + ) — m(D°)\ < 64 MeV in the region p T (K~n + n + ) > 12 GeV and 
\r](K-n + n + )\ > 1.3, or (b) |m(iT~7r + ) - m{D°)\ < 40 MeV elsewhere. Here m(D°) is the 
world average value for the D° mass • This last selection cut is divided into two different 
kinematic regions due to the changing D° mass resolution. The D* + fi~ candidate must 
have an invariant mass in the range 2.5-5.4 GeV. The upper invariant mass cut matches the 
physical requirement of not exceeding the mass of the 5-mesons. 

Because of the kinematics of the D* + decay, the prompt pion takes only a small fraction 
of the energy. The D* + signal is therefore studied as a function of the mass difference Am 
between the D* + and D° candidates. Real D* + mesons are expected to form a peak in Am 
around 145.4 MeV, while the combinatorial background gives a rising distribution, starting 
at the pion mass. The combinatorial background is made of fake D* + yT candidates, created 
from combinations of tracks which pass the selection cuts, but do not come from a D* + 
signal. Figure [I] (a) shows a clear signal in the distribution of Am for the reconstructed 
opposite charge D*fi pairs. The dashed histogram shows the corresponding Am distribution 
for the same charge combinations D* /r* 3 , showing a very small excess around 145.4 MeV, 
whose origin is described in Section [6j 

The opposite charged signal distribution is fitted using a modified Gaussian (G mod ), 
which provides a good description of the tails of the signal distribution. The modified 
Gaussian has the form: 

G mod (a;) oc exp[-0.5 ■ x 1+ ^] (2) 

where x = \(Am — Amo)/a\ and Amo and a, free parameters in the fit, are the mean and 
width of the Am peak. 
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(a) (b) 

Figure 1: (a) Distribution of the mass difference Am for D* /i combinations of opposite charge (points) and 
same charge (dashed line). The solid line shows the result of the fit described in the text, (b) Distribution of 
the opposite charge D*fj, invariant mass, for mass combinations within ±3er of the Am peak, without 
applying the invariant mass cut described in the text. The measured distribution is compared with the 
MC simulation, including the contribution of different sources of signal. The hashed bands show the MC 
statistical uncertainty. 



The combinatorial background is fitted with a power function multiplied by an exponen- 
tial function: 

B(Am) oc (Am - m^e"^"^ (3) 

where a and (3 are free fit parameters, and is the charged pion mass. 

The fitted yield is 4516 ± 100 events, with a fitted Am = 145.463 ± 0.015 MeV, to 
be compared with the world average value 145.421 ± 0.010 MeV [26j, and a fitted o = 
0.49 ± 0.03 MeV. The uncertainties on the fitted Amo and a values are statistical only. 

Figure [l] (b) shows the D* + fi~ invariant mass distribution selected in a region of 3cr 
around the Am peak, without applying any D* + fi~ invariant mass cut. The measured 
distribution is compared with the MC bb+cc simulation described in Section |4j which takes 
into account the contribution of different physical sources to the D* + fi~ signal, as discussed 
in more detail in Section [6] The MC simulation is separately normalised to the number of 
signal and background events in data. The selection on m(D* + fi~) has full efficiency for the 
signal, while rejecting part of the combinatorial background and physical processes other 
than a single 6-hadron decay. 

In order to evaluate differential cross sections, the sample is divided into six p^(D* + 
bins and five \r)(D* + \ bins. The Am distribution in each bin is fitted independently using 
the same fitting procedure as for the total sample. The number of candidates in each bin is 
reported in Table [TJ together with its statistical uncertainty from the fit. 

6. D* + fi~ sample composition 

Various processes contribute to the _D* + /i~ data sample: 
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N(D* + /i-) 


9-12 GeV 


334 ± 33 


12-15 GeV 


1211 ±56 


15-20 GeV 


1527 ±55 


20-30 GeV 


1049 ± 42 


30-45 GeV 


310 ±21 


45-80 GeV 


76 ±10 



\V(D* + H~)\ 




0.0-0.5 


1330 ± 47 


0.5-1.0 


1207 ±47 


1.0-1.5 


919 ±48 


1.5-2.0 


890 ± 60 


2.0-2.5 


317 ±37 



Table 1: Fitted number of opposite charge D* /i pairs for different pt and bins. 



• Direct semileptonic decay: b — > D* fi X; this is the signal contribution used for this 
measurement. 

• Decays of two c-hadrons, one of them decaying semileptonically: c — > D* + X, c — > 

/JL-X'. 

• Direct semileptonic r decay: b — > D* + r~ X, r~ — > ji^u^Uri^). 

• Decays of 6-hadrons with two c-hadrons in the final state, one of them decaying semilep- 
tonically: b ->• D* + DX, D ->• fi'X'. 

• Decays of two 6-hadrons, one of them decaying semileptonically: b — > D* + X, b — > \iX' . 
This source contributes to opposite-sign and same-sign charge combinations, depending 
on the direct or indirect semileptonic decay relative branching ratio and on the neutral 
6-meson oscillation rate. This explains the small excess observed in Figure [l] (a) in the 
peak region of the same sign charge Am distribution. 

• A D* + meson accompanied by a fake muon, contributing to both opposite-sign and 
same-sign charge combinations. The contribution from combinations with misidenti- 
fied muon charge is negligible. 

For the purposes of this measurement, only the direct semileptonic component is of 
interest. Therefore it is necessary to evaluate the fraction of the reconstructed D* + 
sample that actually originates from direct semileptonic b decays. This is estimated from 
the MC simulation. The most significant D* + ^ contributions are listed in Table |2j together 
with the MC statistical uncertainty. 

The fractions from single b semileptonic decays evaluated in the various px and \q\ 
bins of the D* + pair, are reported in Table [3j together with the MC statistical uncertainty 
of the calculations. These values are used for the differential cross section measurements. 



7. Reconstruction and muon trigger efficiency 

The overall efficiency e for Hf, — > D* + fi~X decays to enter the D* + sample, which 
includes the reconstruction, muon trigger and selection efficiencies, is evaluated as a product 
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Source 


Fraction (%) 


b -> D* + n~X 


93.2 ±0.3 


c ->■ c -> /i-X' 


3.8 ±0.2 


6 -> D* + t-X, T~ -> /U-X' 


1.5 ±0.1 


6 D* + DX, D ->■ 


0.9 ±0.1 


others 


0.6 ±0.1 



Table 2: Different sources contributing to the D* /i sample. The uncertainties are due to MC statistics. 





fb (%) 


9-12 GeV 


90.8 ± 1.2 


12-15 GeV 


92.7 ±0.5 


15-20 GeV 


93.8 ±0.4 


20-30 GeV 


93.2 ±0.5 


30-45 GeV 


93.8 ±0.9 


45-80 GeV 


93.1 ± 1.9 



\v(D* + »-)\ 


ft (%) 


0.0-0.5 


93.0 ±0.5 


0.5-1.0 


92.6 ±0.5 


1.0-1.5 


93.4 ±0.6 


1.5-2.0 


93.5 ±0.6 


2.0-2.5 


94.6 ±0.9 



Table 3: Fractions of single b semileptonic decays in different pr(D* /z ) and \rj(D* jji )| bins. The 
uncertainties are due to MC statistics. 



of three different components, in order to combine MC and data-driven efficiency calcula- 
tions. Since the only requirement is the single b detection efficiency, the bb MC sample is 
used. The components are defined as: 



e trigg 



^selection 



iV(true D* + \i with /i and tracks reconstructed) 
6rec ° ~ iV(true D* + n~) 

iV(true D* + with fi and tracks reconstructed, \i matched to trigger) 



(4) 



iV(true D* + fi~ with ji and tracks reconstructed) 
iV(true D* + with /i and tracks rec, \i matched to trigger, D* + ^ selection 
iV(true D* + fi~ with \x and tracks rec, \x matched to trigger) 



(5) 



(6) 

where the number of true D* + fi pairs is calculated within the fiducial kinematic region 
Pt(D* + ) > 4.5 GeV, p T (/0 > 6 GeV, \r](D* + )\ < 2.5 and \rj(v~)\ < 2.4. Events where 
the D* + is inside the fiducial region, but its decay products are not fully reconstructed, 
contribute to e reco - 

Both e roco and ejection are taken from MC simulation. However etri gg er, which is the 
fraction of the reconstructed muons that actually satisfied the trigger, is measured directly 
from data using J/ip — > fi + n~ samples [42] . These efficiencies are evaluated for the same 
data-taking periods used in this measurement. 

The overall efficiency e is given by: 



e — e rcco(MC)et r ig g er(data)e se i ec tion(MC) 
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(7) 



The different efficiency components, together with the related statistical uncertainties, 
are determined as e reco = (48.3 ±0.4)%, e tr igg er = (81.9 ±0.4)% and e se i cction = (79.1 ±0.5)%. 
The overall efficiency is (31.3±0.4)%, and the values obtained in p^(D* + and \r}(D* + fi~)\ 
bins are reported in Table |4} A complete description of the systematic uncertainties follows 
in Section [U 





e(%) 


9-12 GeV 


21.2 ±0.9 


12-15 GeV 


26.7 ±0.6 


15-20 GeV 


32.1 ±0.6 


20-30 GeV 


38.8 ±0.9 


30-45 GeV 


45.2 ± 1.7 


45-80 GeV 


52 ±4 





e(%) 


0.0-0.5 


37.5 ±0.7 


0.5-1.0 


37.2 ±0.8 


1.0-1.5 


29.9 ±0.8 


1.5-2.0 


26.1 ±0.8 


2.0-2.5 


16.1 ±0.9 



Table 4: Overall efficiency e for different pt(D* + [i ) and \r](D* + /i )| bins. 



8. Systematic uncertainties 

The uncertainty in the cross section due to each systematic variation is evaluated by 
repeating the entire analysis procedure and finding the change in the cross section value. The 
same strategy is adopted to evaluate bin-by-bin systematic uncertainties for the differential 
cross section measurements. The following sources are considered: 

• Uncertainty of the yields from the fits, obtained by varying the fitting procedure in 
the following ways: 

— reducing the high end of the Am range used for the D* + fi~ signal fit by 4MeV, 
from 165 MeV to 161 MeV; 

— changing the background parameterisation function to be oc 1 — exp(— a (Am — 
nT>ir) )i where a and are free fit parameters, which provides a P(x 2 ) for the fit 
similar to that with the default background parameterisation. 

• Uncertainty of the sample composition estimate: the ft, measurement depends on 
the b/c cross section ratio used in the MC sample. The ratio of the beauty and 
charm contributions to the inclusive D* + production, estimated using the life-time 
information, has been found to be in agreement with the ratio in Pythia, within 
experimental uncertainties. To cover the uncertainties, the MC b/c ratio is varied 
between 50% and 200% of its nominal value. 

• Uncertainties of the muon trigger efficiencies are estimated from J /tp — > ji + pT studies 
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Uncertainty on tracking and muon reconstruction efficiency: the uncertainty on ID 
tracking efficiency is dominated by the detector material description used in MC sim- 
ulations. This uncertainty is evaluated in studies of minimum bias events [28]. The 
muon reconstruction uncertainty is evaluated on Z — > fi + fi~ data samples [12] • This 
systematic uncertainty is dominated by the ID tracking uncertainty. 

Model dependence of the reconstruction efficiency: the efficiency calculation could be 
affected by differences between the pt(-D* + aO and r](D* + [i~) spectra in data and MC 
simulation. To estimate the systematic uncertainty, the MC distribution is varied, 
while preserving consistency with the observed data distribution, and the resulting 
change in efficiency is computed after each variation. 

Uncertainty due to differences in the fit of the D° and 6-hadron vertices between data 
and MC simulation: to estimate the systematic uncertainty, the MC P{x 2 ) distribution 
is varied, while preserving consistency with the observed data distribution, and the 
resulting change in efficiency is computed after each variation. 

Uncertainty of the difference in D° mass resolution between data and MC simulation: 
the efficiency calculation is corrected to account the difference between D° mass reso- 
lution in data and MC simulation. To estimate the systematic uncertainty, the error 
on the data-to-MC ratio of D° mass widths is propagated to the efficiency. 

NLO prediction uncertainty: since the NLO predictions are also used as an active 
part of the analysis for unfolding (Section [9]) and acceptance corrections (Section 10) 



the theoretical uncertainties and the use of different predictions introduce additional 
systematic uncertainties to the experimental measurements. These are evaluated by 
repeating the entire analysis, introducing different theoretical uncertainties (Section [4]) 
to the default central prediction (Powheg+Pythia), and using a different theoret- 
ical prediction (POWHEG+HERWIG and MC@NLO): positive and negative differences 
obtained with respect to using the central prediction are separately summed in quadra- 
ture. The use of the predictions matched with Herwig produces visible asymmetries 



in the uncertainties of the acceptance corrections (Section 10). 
Uncertainty of the luminosity measurement (±3.4%) [33 



Relative uncertainty on the branching fractions of the different decay chains, obtained 
from the world averages [26]: & ^ D* + (±7 %), D* + D°ir+ (±0.7 %), D° -> 
K-n+(±l.S%). 



In Sections [9] and 10 , tables are shown with these uncertainties quoted after each step of 
the analysis. 
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9. Differential cross sections for H b — >• D* + n X production 

Differential cross sections for Hf, — > D* + fi~X production as a function of the and \rj\ 
of the D* + fi~ pairs are evaluated by using Equation [l] and dividing by the bin width. The 
results are shown in Table [5] 



[GeV] 


da(H b ^D*+^-X) 
dp T (D*+^-) 

[nb/GeV] 


9-12 


2.78 ± 0.29±g;|g 


12-15 


8.2 ± 0.4±g;| 


15-20 


5.2±0.2l°j 


20-30 


1.47 ± 0.06±g;i| 


30-45 


0.250 ± 0.018±g;g|| 


45-80 


0.0229 ± 0.0030±g:g^| 



\v(D* + t*-)\ 




d[7?(D*+A t -)| 

[nb/unit of \rj\] 


0.0-0.5 


38.4 ± 1-5±|^ 


0.5-1.0 


34.9±1.4j|;J 


1.0-1.5 


33.5 ± 1.8±|;f 


1.5-2.0 


37.2 ±2.6^ 


2.0-2.5 


21.7 =b 2.6±|;[ 



Tabic 5: Differential cross sections for Hb — > D* + fj,~ X production as a function of pt and \rj\ of the 
D* + pair, in the fiducial kinematical region pt(D* + ) > 4.5 GeV,pr(/-t _ ) > 6 GeV, \r](D* + )\ < 2.5 and 
|?7(/i~)| < 2.4. The statistical and total systematic uncertainties are shown for each cross section. 



To extract differential cross sections as a function of the pr and \rj\ of the 6-hadron, it 
is necessary to correct the observed pt(D* + fi~) and \i](D* + fi~)\ distributions using Monte 
Carlo simulations, in order to take into account the kinematics of the missing particles from 
the decay Hb — > D* + ^rX. This procedure is known as unfolding [4"5] - [4"8] . The unfolding 
approach used in this paper is based on the iterative method described in Ref. [49J, containing 
elements of Bayesian statistics. 

The element of the response matrix F for a 6-hadron in a Pt/M(-H&) bin j to decay 
into a D* + in p^/\r]\(D* + fi~) bin % can be interpreted as a conditional probability 

= P(D* + ^ in bin i\H b in bin j). (8) 

Given an initial set of probabilities pi for 6-hadrons to be found in bin i, using Bayes' 
theorem one can obtain the expected number of 6-hadrons in bin i, given a measured D* + 
distribution: 

N** b = p ( H b m bin i\D*n in bin j)Nf* + »~ = £ f j ^ )Nj >m+t '~ . (9) 

A NLO Monte Carlo sample generated with Powheg+Pythia is used to create the default 
response matrix F and the initial prior probabilities p. The procedure is repeated with 
different MC generators, in order to evaluate systematic uncertainties. 

The procedure can be iterated, taking as new prior probabilities the solutions of the pre- 
vious step, i.e. pi = N^/N^. After a large number of iterations, the procedure converges 
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on the results obtained with a direct inversion of the response matrix F 

AW 
3=1 

This method is known to be sensitive to statistical fluctuations [45J, but this effect can be 
mitigated in the Bayesian method by truncating the procedure after a few iterations. 

The number of iterations was therefore optimised in Monte Carlo simulations with test 
measurements, comparing the values obtained after each iteration to the values expected 
from the MC-generated information, using a \ 2 test. Two iterations are the optimal solution 
in this case, providing compatible results even when the response matrix F and the prior 
probabilities p are generated using different theoretical distributions. 

The inversion method and the Bayesian method with a different number of iterations 
were employed as a check. Within the systematic uncertainties, all the results were found 
to be in agreement with the chosen default procedure. 

A bias could occur in this procedure due to the possible mismodelling of the Hb decays 
(e.g. D** decays contributing to the missing particles in the final state) in the simulation. It 
was verified with the simulation that the relevant D* + kinematic variables have a small 
dependence on the specific 6-hadron decay, and that a mismodelling of the D** branching 
ratios does not produce a significant effect. This is expected since the dominant D* + jjT 
contribution arises from direct B° decays without an intermediate D** . 

Once the Hb distribution is obtained, the differential Hb — > D* + fi~X cross sections 
are determined as a function of pt and \rj\ of the 6-hadron, inside the kinematic region 
p T (D* + ) > 4.5 GeV, p T (/0 > 6 GeV, \rj(D* + )\ < 2.5 and \r)(jT)\ < 2.4. 

Figure [2] shows the measured differential cross sections, with comparisons to the NLO 
theoretical predictions. The Powheg+Pythia shaded band refers to the total theoretical 
uncertainty of the prediction. The differential cross section values are reported in Table |6j 
together with the statistical and total systematic uncertainties. The individual contributions 
to the systematic uncertainties are listed in Tables [7] and [8] The comparison with data shows 
that NLO calculations underestimate the cross section, although the difference is within the 
combined experimental and theoretical uncertainties. 

The integrated Hb — > D* + X cross section, inside the kinematic region p^{D* + ) > 
4.5 GeV, Pt{h~) > 6 GeV, \r](D* + )\ < 2.5 and \r)(p-)\ < 2.4, is: 

a(pp H b X' D* + fi-X) = 78.7 ± 2.0 (stat.) ± 7.3 (syst.) ± 1.2(^) ± 2.7(^f ) nb 

The integrated Powheg+Pythia prediction, with its theoretical uncertainty, is: 

a(pp -> H h X' D* + /j-X) = 53i^(scale) ±l(m b ) ±3 (PDF) ^(hadr.) nb 

The corresponding Powheg+Herwig prediction is 51 nb, while mc@nlo predicts 56 nb, 
with similar theoretical uncertainties to the Powheg+Pythia prediction. 
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Figure 2: Differential cross section for Hb — > D" (J>~ X production as a function of (a) pt and (b) \r]\ of 
the &-hadron, in the fiducial kinematical region p^(D* + ) > 4.5 GeV,pr(M~) > 6 GeV, \rj(D* + )\ < 2.5 and 
|?7(^ _ )| < 2.4. The measurement is compared with the theoretical predictions, as described in the text. The 
inner error bars of the data points are statistical uncertainties, the outer are statistical+total systematic 
uncertainties. 



10. Differential cross sections for b-hadron production 

The fe-hadron differential cross sections can be derived from the H b — > D* + fi~X differ- 
ential cross sections by taking into account the branching ratio £${b — > D* + ^rX) and the 
necessary decay acceptance corrections. These are evaluated using a Powheg+Pythia 
simulation in two steps: 

• Identification of the H b kinematic region selected by the D* + and /i~ kinematic cuts. 
This indicates that only 6-hadrons with p^(H b ) > 9 GeV and \r](H b )\ < 2.5 pass the 
D* + and \T kinematic cuts. 

• Evaluation of a bin-by-bin p^- and |?7|-decay acceptance a in the H b allowed kinematic 
region, defined as 

number of H b (— > D* + fi~) passing the D* and fi kinematic cuts 

a = t — : (11) 

number of #&(—»■ D* passing the H b kinematic cuts 

The results are shown in Table [9] for the Powheg+Pythia central prediction. Section [8] 
describes how the NLO theoretical uncertainties are propagated to this measurement. 

The 6-hadron differential cross sections as a function of pt and rj, inside the kinematic 
region p T (H b ) > 9 GeV and \r](H b )\ < 2.5, can then be calculated according to the formula: 

da(H b X) = 1 da(pp^H b X' ^D*+(i-X) 

dp T (r]) " a PTM ^(b -> D* + pi~X) dp T (r}) 
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PT{H b ) 

[GeV] 


&a{H b -+D*+ X) 


da(H b X) 


[nb/GeV] 


[nb/GeV] 


9-12 


0.73 ±0.12lg;?? 


(5.8 ±0.9l?;{j) • 10 3 


12-15 


4.65 ± 0.27^lo 


(2.37 ±0.14l8;l) • 10 3 


15-20 


5.48 ±0.191811 


(9.1 ± 0.31^) -10 2 


20-30 


2.46 ± 0.08l8;24 


212±7 +2 | 


30-45 


0.530 ± 0.025l8:8e2 


31.3 d= 1-5±1;| 


45-80 


0.055 ± 0.005l8:88e 


2.78 ± 0.2518:1 



\v(H b )\ 


dcr(H b ^D*+/i-X) 


da(H b X) 


d\v{H b )\ 

[nb/unit of \r]\] 


d\v{H b )\ 

[/xb/unit of \r)\] 


0.0-0.5 


38.0 ± 1.51 3 ; 3 


14.3 ± 0.61^ 


0.5-1.0 


35.0 ± 1.5±|;1 


13.4 ± 0.61^ 


1.0-1.5 


32.9 ± 1-9±1:? 


13.1 ± 0.711;^ 


1.5-2.0 


37.5 ± 2.7ttl 


i5.8 ± 1.1+4:! 


2.0-2.5 


22.3 ± 2.8+1:1 


13.3 ± 1.6 +2 J 



Table 6: Differential cross sections for Hi, —> D* + [i X and HhX production as a function of px and |^| 
of the fr-hadron, in the fiducial kinematical regions pt (D* + ) > 4.5 GeV.pr&O > 6 GeV, \rj(D* + )\ < 2.5 
and \r](iJ,~)\ < 2.4, and pT(Hb) > 9 GeV, \r](Hb)\ < 2.5 respectively. The statistical and total systematic 
uncertainties are shown for each cross section. 

Figure [3] shows the 6-hadron differential cross section measurements compared with the- 
oretical predictions. The shaded band is the overall theoretical uncertainty of the central 
Powheg±Pythia prediction. Since the acceptance correction factors have a dependence 
onpr and \r]\, as shown in Table |9j the shapes of the 6-hadron differential cross sections are 
different to the H b — > D* + fi~X differential cross sections shown in Figure [2j The systematic 
uncertainties are those from the a(H b — > D* + fi~ X) measurement described in Section |9j 
with the addition of the uncertainty of the branching ratio &(b — > D* + fi~X) and the un- 
certainties of the decay acceptance correction. The 6-hadron differential cross section values 
are reported in Table [6j together with the statistical and total systematic uncertainties, 
while the individual contributions to the systematic uncertainty are reported in Tables [7] 
and [8j The combined unfolding and acceptance uncertainties are calculated taking their 
correlations into account. 

The comparison with data shows that NLO calculations underestimate the cross section, 
although the difference is within the combined experimental and theoretical uncertainties. 
The 6-hadron integrated cross section for pi{H b ) > 9 GeV and \rj(H b )\ < 2.5 is measured as: 

a(pp^ H b X) = 32.7 ±0.8 (stat.) ± 3.1 (syst.) +2 ,;J(a) ± 2.3(^) ± l.l(jSf) fib 

The integrated Powheg±Pythia prediction, with its theoretical uncertainty, is: 
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Pt bin (GeV) 


9-12 


12-15 


15-20 


20-30 


30-45 


45-80 


data statistics 


±15.8 


±5.9 


±3.4 


±3.1 


±4.7 


±9.0 


a(H b — > D* + fi X) and a(H b ) relative systematic error (%) 


D*fi fit 


±3.5 


±1.8 


±1.0 


±1.4 


±1.7 


±2.0 


fb 


+2.5 
-3.8 


+2.3 
-3.5 


+1.8 
-2.8 


+1.6 
-2.5 


+1.4 
-2.2 


+1.8 
-2.9 


H trigger 


+1.3 
-1.2 


+1.3 
-1.3 


+ 1.7 
-1.6 


+2.2 
-2.0 


+2.5 
-2.2 


+2.7 
-2.5 


tracking + fi reconstruction 


+9.1 
-8.2 


+9.0 
-8.1 


+8.9 
-8.0 


+8.7 
-7.9 


+8.5 
-7.7 


+8.3 
-7.5 


MC Pt/v reweight 


+0.2 
-1.3 


+0.2 
-1.2 


+0.4 
-1.1 


+0.5 
-1.1 


+0.4 
-1.0 


+0.2 
-0.8 


D° and H b vertices fit 


±2.0 


±2.0 


±2.0 


±2.0 


±2.0 


±2.0 


D° mass correction 


+0.8 
-1.0 


+0.8 
-1.0 


+0.8 
-1.0 


+0.8 
-1.0 


+0.8 
-1.0 


+0.8 
-1.0 


luminosity 


±3.4 


^{D* + -> D°7i + ) 


±0.7 


^{D° K--K+) 


±1.3 


a(H b — > D* + fi X) relative systematic error (%) 




unfolding 


+6.6 
-10.0 


+2.3 
-3.8 


+1.7 
-1.5 


+2.3 
-1.3 


+3.2 
-6.7 


+9.1 
-3.5 


<j(H b ) relative systematic error in (%) 


3S{b D* + fi~X) 


±7 


unfolding © acceptance 


+3.4 
-11.3 


+1.6 
-6.4 


+2.3 
-4.0 


+0.7 
-2.5 


+1.7 
-4.4 


+6.0 
-1.0 


total syst a(H b — > D* + \T X) 


+ 12.9 
-14.7 


+10.7 
-10.8 


+10.3 
-9.9 


+ 10.4 
-9.8 


+10.6 
-11.7 


+13.6 
-10.3 


total syst a(H b ) 


+ 13.4 
-17.1 


+12.6 
-14.0 


+12.5 
-12.5 


+ 12.3 
-12.1 


+12.3 
-12.5 


+13.5 
-11.8 



Table 7: Hb — > D* + /i X and cross section relative uncertainties as a function of pT(Hb), listed as 
percentages (%). 



a(pp ->• H b X) = 22.2+^(scale) t?;J(m 6 ) tfi(PDF) -Ls(hadr.) fib 

The corresponding Powheg+Herwig prediction is 18.6 fib, while mc@nlo predicts 
19.2 /ib, with similar theoretical uncertainties to the Powheg+Pythia prediction. 

11. Discussion 



Section 10 discusses the measurement of the 6-hadron production cross section for p^{Hb) > 
9 GeV and \r}{Hi)\ < 2.5. In order to compare this result with other LHC measurements, we 
extrapolate this measurement to the full kinematic phase space, extending to regions out- 
side the ATLAS coverage, using the NLO MC theoretical predictions. The multiplicative 
extrapolation factor is defined as the ratio of the total number of generated 6-hadrons to the 
number of 6-hadrons generated withpT(^) > 9 GeV and \r](Hb)\ < 2.5, and is estimated to 
be ll.Oj^g. The resulting total 6-hadron cross section is: 

cr(pp -»■ H b X) total = 360 ± 9 (stat.) ± 34 (syst.) ± 25(^) ± 12 (if) ^(accept. © extrap.) ^b 
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-6.5 
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+8.5 
-7.7 
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+16.2 
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+8.6 
-8.0 


MC Pt/v reweight 


+1.5 
-0.1 


+1.2 
-0.1 


+1.4 
-0.1 
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-0.2 
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D° and Hb vertices fit 
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-1.0 


+0.8 
-1.0 
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-1.0 
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±3.4 
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-1.5 


+1.4 
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cr(Hb) relative systematic error (%) 


38 {b -> D* + fi~X) 


±7 


unfolding © acceptance 


+5.1 
-15.0 


+7.3 
-16.2 


+10.7 
-19.1 


+4.8 
-24.6 


+4.3 
-29.6 


+6.4 
-17.1 


total syst a(H b D* + fi~X) 


+8.8 
-8.5 


+9.0 
-9.0 


+10.0 
-9.4 


+12.5 
-11.4 


+17.1 
-14.5 


+10.0 
-9.8 


total syst cr(Hb) 


+12.2 
-18.5 


+13.4 
-19.7 


+16.1 
-22.3 


+15.0 
-27.9 


+18.8 
-33.5 


+13.8 
-20.9 



Table 8: Hi, — > D* + /i X and Hb cross section relative uncertainties as a function of \r](Hi,)\, listed as 
percentages (%). The last column refers to the integrated cross sections. 

where the combined acceptance and extrapolation uncertainty is calculated taking their 
correlations into account. 

This value can be compared with the inclusive bb cross section measurements by LHCb 
a(pp —¥ bbX) = 284 ±20 (stat.) ±49 (syst.) fib, evaluated in the kinematic region 2 < r\ < 6 
using decays to D^iTvX final states [IE], and a{pp —> bbX) = 288 ±4 (stat.) ±48 (syst.) fib, 
evaluated using J /ipX final states in the kinematic region 2.0 < y < 4.5 [17]. Extrapolations 
outside the LHCb sensitivity region are done using different theoretical models, without 
including additional uncertainties. Also ALICE measured the inclusive bb cross section 
in pp collisions, using decays to J/ipX final states in the kinematic region \y\ < 0.9 and 
Pt > 1-3 GeV [21]. After extrapolation to the full phase space, they obtain a{pp — > bbX) = 
244 ±64 (stat.)^S (syst.)^(extr.) fib. 
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PT(H b ) 


a 


9-12 GeV 


0.005 


12-15 GeV 


0.071 


15-20 GeV 
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30-45 GeV 


0.614 
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1.5-2.0 


0.086 


2.0-2.5 
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Table 9: Decay acceptance a as a function of pT{H b ) and \rj(H b )\ for the Powheg+Pythia prediction. 
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Figure 3: Differential cross section for H b production as a function of (a) pt and (b) \rj\ of the 6-hadron, 
in the fiducial kinematical region p^{H b ) > 9 GeV, |ry(ffj,)| < 2.5. The measurement is compared with the 
theoretical predictions, as described in the text. The inner error bars of the data points are statistical 
uncertainties, the outer are statistical+total systematic uncertainties. 



12. Conclusions 

The production of 6-hadrons (Hb) at the LHC is measured with the ATLAS detector in 
proton-proton collisions at ^/s = 7 TeV, using 3.3 pb -1 of integrated luminosity from the 
2010 run. A 6-hadron enriched sample was obtained by combining oppositely charged D* 
mesons and muons, in events triggered by a muon with pp exceeding 6 GeV. 

Differential cross sections as functions of pt and \rj\ are produced for both H b and 
H b —7- D* + fi~X production. These measurements are found to be higher than the NLO 
QCD predictions, but consistent within the experimental and theoretical uncertainties. The 
integrated 6-hadron cross section for p T (H b ) > 9 GeV and \r)(H b )\ < 2.5 is measured as 

a(pp-> H b X) = 32.7 ±0.8 (stat.) ±3.1 (syst.)^;J(a) ± 2.3(30) ± l.l(jSf) fib. 
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